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Abstract 

We present an equation of state for elastic matter which allows for purely longitudinal elastic 
waves in all propagation directions, not just principal directions. The speed of these waves is 
equal to the speed of light whereas the transversal type speeds are also very high, comparable 
to but always strictly less than that of light. Clearly such an equation of state does not give a 
reasonable matter description for the crust of a neutron star, but it does provide a nice causal 
toy model for an extremely rigid phase in a neutron star core, should such a phase exist. Another 
reason for focusing on this particular equation of state is simply that it leads to a very simple 
recipe for finding stationary rigid motion exact solutions to the Einstein equations. In fact, we 
show that a very large class of stationary spacetimes with constant Ricci scalar can be interpreted 
as rigid motion solutions with this matter source. We use the recipe to derive a static spherically 
symmetric exact solution with constant energy density, regular centre and finite radius, having 
a nontrivial parameter that can be varied to yield a mass-radius curve from which stability can 
be read off. It turns out that the solution is stable down to a tenuity R/M slightly less than 3. 
The result of this static approach to stability is confirmed by a numerical determination of the 
fundamental radial oscillation mode frequency. We also present another solution with outwards 
decreasing energy density. Unfortunately, this solution only has a trivial scaling parameter and 
is found to be unstable. 

PACS: 04.40.Dg, 97.10.Cv, 97.60.Jd 



1 Introduction 

There are few exact solutions to the Einstein equations that can be considered stellar models. The most 
well-known solution is the interior Schwarzschild solution which can be thought as a limiting case for perfect 
fluids since it satisfies the Buchdahl inequality R/M < 9/4 as an equality when the central pressure is 
taken to infinity. However, the constant energy density of this solution makes it unphysical, at least when 
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interpreting it as an adiabatic perfect fluid since the speed of sound is then infinite. Ahhough there are 
other known perfect fluid solutions, only a few of them satisfy the basic criteria of physicality, cf. for 
a comprehensive classification of known exact solutions. Several workers have also considered equilibrium 
models with pressure anisotropy, at least since the paper 0. In many cases an "exact solution" is obtained 
through some ansatz which does not correspond to a particular matter description. Elastic matter, on 
the other hand, gives rise to anisotropic pressures in a natural way and should moreover be relevant for 
modelling of e.g. neutron stars since they are believed to have solid crusts and, perhaps, solid cores. The 
general relativistic theory of elasticity has never been a very hot field of research, but over the years it 
has generated quite a number of papers, most of which concern the theoretical formulation of the theory. 
Despite of this, the times when the theory - in its full nonlinear regime - has found its way into astrophysical 
applications are scarce. 

In a recent paper by the present authors[3] (hereinafter paper I), we reconsidered the subject and applied 
it to static spherically symmetric (SSS) configurations. In particular we constructed elastic stellar models 
numerically for a specific equation of state. In a subsequent paperQ (paper II) we showed how to analyse the 
radial oscillations of elastic stars and found, in particular, that the models constructed in paper I are stable 
up to the first maximum of the mass-radius curves, as one would expect. Although one will probably always 
be forced to resort to numerical methods when one wants to model a real star, it would nevertheless be useful 
to have a simple exact solution that is not too unphysical. Such a solution would serve as a convenient tool 
for investigating qualitative features and for testing numerical codes. In this context it could be useful to 
list some criteria that an exact SSS solution should satisfy in order to be of astrophysical interest. In our 
humble opinion, this list should minimally include 

• The centre should be regular (elementary flat), implying pressure isotropy there. 

• The energy density p should be positive at the centre and everywhere non-negative. 

• The isotropic central pressure pc should be positive. 

• The surface of the model, where the radial pressure Pr by necessity is zero, should be at a finite 
Schwarzschild radius r = R. Usually, although perhaps not necessarily, there are no other zeros of pr 
at smaller radii, implying that pr > everywhere inside the surface of the star. 

• The three elastic wave modes associated with any propagation direction should have speeds v obeying 
< < 1 in geometrical units. In other words the equation of state, at least in the range it is used, 
should be microstable (u^ > 0) and causal (u^ < 1). 

• The model should be stable against radial perturbations. 

Exact solutions with elastic matter sources have previously been constructed by Magli and Kijowski[n]. 
These solutions do not obey all of the first four criteria however, while to our knowledge it is undetermined 
whether they obey the last two. On the other hand it should be stressed that a solution may fulfill the 
criteria only when matched to some other matter region instead of used as a stellar model in its own right. 
In this paper we will present a family of solutions having a branch on which all of the criteria are satisfied. 
The contents and structure of the present paper are as follows; 

In section 121 we introduce and discuss the properties of a particular equation of state, which we refer to 
as the stiff ultrarigid equation of state (SUREOS). The most important features of this equation of state are 
listed below. 

• The SUREOS can be split into a conformally invariant part and a vacuum energy density leading 
to the relation p = pi + P2 + P3 + 4i? for the energy density and principal pressures which in turn 
implies that the spacetime Ricci scalar is constant and equal to AnB, with k = Stt in geometrical units. 
In this sense the SUREOS is similar to the MIT bag (perfect fluid) equation of state p = 3p -I- AB 
describing non-interacting strange quark matter in the limit of zero quark massjHI. Furthermore, the 
energy density and principal pressures are required to satisfy the inequality p > p^ + 2B for all p. 
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• The SUREOS is stijf in the sense that the speed of sound for longitudinal elastic waves propagating in 
principal directions is marginally causal, i.e. equal to the speed of light. In fact the SUREOS allows for 
one purely longitudinal wave mode with marginally causal sound speed in any propagation direction. 
In this sense the SUREOS is similar to the stiff perfect fluid equation of state p = p + po proposed by 
Zel'dovichQ. 

• The SUREOS is ultrarigid in the sense that it has a very high shear modulus, of the same order as the 
energy density. This results in very high propagation speeds for the remaining two elastic wave modes 
that are purely transversal in principal directions. In an unsheared state these speeds are 1/^/2 times 
the speed of light. In general they are always real, positive and strictly less than the speed of light. 

In section 01 we give a recipe for obtaining stationary rigid motion solutions to the Einstein equations 
with SUREOS elastic matter source. The recipe simply amounts to making the spacetime Ricci scalar take 
a constant value to be identified with AkB and, in addition, to ensure that the inequalities p > + 2B are 
satisfied by the eigenvalues of the resulting stress-energy tensor. As soon as this has been achieved one has 
an exact solution to the Einstein equations with SUREOS elastic matter source, with respect to a material 
space metric that can be explicitly constructed according to a simple algebraic formula. Clearly one can 
take the view that this method of letting the solution determine the material space metric, rather than the 
other way around, is cheating. Indeed, should one want to specify the material space metric by e.g. taking 
it to be flat, then the method should be of little use since the resulting "backwards" constructed metric 
will in general be curved. Although a flat metric is the correct choice when one wants to describe a solid 
region without any lattice dislocations (c/. section 3 in paper I), such a desciption will only be appropriate 
throughout a rather small solid region of e.g. a neutron star. In particular, a neutron star crust, having a 
thickness of the order 1 km, is bound to contain many dislocations. Unless one wants to go into a detailed 
description of the dislocations, their average effect over distances containing many of them can instead be 
considered to be encoded into an effective, generally curved, material space metric. In the case of neutron 
star crusts this effective metric is in principle "created" at the moment when the outer region of a new 
born neutron star cools down to the temperature at which solidification occurs, or at the moment when an 
older/colder neutron star settles down after a star quake (rupture of the crust). Analogous remarks apply 
to solid cores in neutron stars, should such cores exist. When constructing an equilibrium stellar model, 
the simplest choice of effective material space metric is the one that results (usually only implicitly) from 
assuming that the star is in an unsheared state, implying pressure isotropy so that the star can alternatively 
be thought of as being made up out of perfect fluid matter. In this case the rigidity of the stellar matter only 
enters through a nonvanishing shear modulus when considering stellar perturbations around the equilibrium 
model. Whereas this approach is fine under many circumstances, it does not allow for studies of how the 
perturbations are affected by pressure anisotropics of the background. To test the effect of anisotropics one 
can use moderately (and hence realistically) anisotropic equilibrium models, such as the ones constructed 
numerically in paper I, or some exact solution which is likely to have more extreme pressure profiles, should 
it be possible to find one that obeys the basic list of criteria set up above. This paper goes to show that this 
indeed is possible. 

In section 0] we apply the general recipe of section |31 to construct a spherically symmetric family of 
solutions for which each member has constant energy density. As a direct consequence of the constant energy 
density and the equation of state, the radial and tangential pressure cannot be simultaneously decreasing 
with increasing radius because of the linear relation p = pr + 2pt + AB. Instead it turns out that the 
radial pressure is monotonically decreasing from the centre to the surface while the tangential pressure is 
monotonically increasing at half the rate. For each fixed value of the equation of state parameter B the 
solution family has one remaining free parameter which can be chosen to be the (isotropic) central pressure, 
going between zero and infinity. For 5 = the equation of state is scale invariant and the central pressure 
is just a trivial scaling parameter. However, for i? > (we see < as unphysical) the freely specifiable 
central pressure results in a nontrivial mass-radius curve, exhibiting a mass maximum. From the conclusions 
drawn in paper II this maximum should be the point where instability sets in. More precisely, since the only 
mass extremum is for the maximum mass solution, the solutions with lower central pressure should be stable 
whereas those with higher central pressure should have exactly one unstable mode of radial oscillations, 
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given that it is known that the models are stable in the limit of zero central pressure. This behaviour is 
exactly what a numerical determination of the first two frequencies of radial oscillations shows. The main 
motivation for this paper is the fact that the solutions on the stable branch of the mass-radius curve satisfy 
all items on the above list of criteria for physicality of exact spherically symmetric solutions. Due to the 
rather extreme equation of state and the nature of the energy density and pressure profiles, we make no 
claims that this solution family provides an accurate description of real stars. However, solid neutron star 
cores having very high shear moduli up to the same order of magnitude as the pressure have actually been 
discussed in the literature, especially in the first half of the 70s (c/. Haensel|Hl for a discussion, historical 
notes and further references). Should such cores exist, despite the fact that they have not been theoretically 
favoured in later years, the SUREOS would be a very relevant toy model equation of state. Moreover, we 
wish to emphasize that in contrast to the interior Schwarzschild models which also have constant energy 
density, the elastic matter models presented here have causal speeds of wave propagation. The reason why 
causality is compatible with constant energy density in the latter case is because the energy density consists 
of two parts; one purely compressional part which is outwards decreasing and one shearing part which is 
outwards increasing due to increase of pressure anisotropy. Clearly then, a constant energy density can be 
obtained if these two effects exactly cancel each other out, which is precisely what happens here. 

In section |S1 we present another spherically symmetric solution with B = for which the energy density 
as well as the radial and tangential pressures all are outwards monotonically decreasing, which in that sense 
makes it a more realistic model of a star. The solution turns out to be unstable, however, which should not 
be too surprising considering that a B = model with finite central pressure can be viewed as a rescaled 
version of a B 7^ model with infinite central pressure. Because of its instability we do not discuss this 
solution in great detail, but it should be interesting to generalise it to i? > since we then expect there to be 
a stable branch of solutions anologous to the stable branch of the constant energy density family. The reason 
why there exists more than one one-parameter solution for a given equation of state and central pressure is 
due to the freedom of choosing the material space metric in more than one way. Clearly it is this freedom 
that makes the simple recipe for generating SUREOS solutions possible. Even more interesting would be to 
use the recipe for finding rigidly rotating solutions, but that is beyond the scope of the present paper. 

The notation and conventions will follow papers I and II. 



2 The stiff ultrarigid equation of state (SUREOS) 

As shown in paper I, given an elastic equation of state of the general type associated with a material space 
metric, 

p = p(ni,n2,n3), (1) 

where the n^'s are the principal linear particle densities as defined in paper I, the speed of a longitudinal 
elastic wave propagating in the direction of a unit principal vector is given by 

Using that the principal pressure is defined through the relation 



we can rewrite eq. Q as 



n2 



- dn,. 
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Setting w^ii to a constant, say C||, it follows that p should satisfy the differential equation 



2 9V _ 2 dp 

Since we are only considering elastic materials that are intrinsically isotropic, the equation of state should 
be invariant under permutations of the n^'s. With this restriction the general solution to eq. (jSJ is 

p = Ci{nin2n^T + C2 [(nma)" + (nsng)" + (rigni)"] + C^{n^ + + <) + Ci, a = 1 + C|f , (6) 

where Ci - C4 are constants. In this paper we focus on a particular equation of state belonging to this 
class, namely the one specified by setting a = 2, Ci = C3 = which we refer to as the stiff ultrarigid 
equation of state (SUREOS) for reasons explained in the introduction and below. Renaming the remaining 
two constants as C2 — A, ^ B we find that the energy density and principal pressures for the SUREOS 
are 

p^A [{nni2f + {n2n:if + {n^mf] + B (7) 
dp 

= "m7^ P = ^ W{n^+i + n^+2) - {n^+inf,+2f] - S, (8) 

where we use the cyclic rule /i + 3 = /i. This leads to the linear relation 

p = pi +P2 +P3 + 4-B. (9) 

In an isotropic state, i.e. when ni = ^2 = 713 = 'n}/'^ leading to Pi = P2 = Ps = P, this relation mimics 
the MIT bag equation of state p — ip + AB with B being the bag constant. This directly suggests that B 
should be positive or at least non-negative, since B < Q would lead to a negative energy density at zero 
isotropic pressure. Now, while the speed of sound for the MIT bag equation of state is times the speed 

of light, the longitudinal speed of sound in principal directions for the SUREOS is by construction equal to 
the speed of light since we have set a = 2. With the longitudinal speed of sound being marginally causal 
in principal directions, it is natural to wonder if it may become acausal in other propagation directions. In 
general, neither longitudinal nor transversal waves stay purely longitudinal/transversal when continuously 
varying the propagation direction away from a principal direction, but the SUREOS in fact allows for a 
purely longitudinally polarized wave mode in all propagation directions. Moreover, the speed of these waves 
is direction independent and hence equal to the speed of light. These statements are proved in the appendix. 
What about the remaining two polarization modes? As opposed to the longitudinal mode their speeds 
generally depend on the polarization direction. In principal directions the modes are purely transversal with 
speeds v^_\_v given by 

^/fi 1/ = (p + Piy)~^ " ^'^ 7^— — — TT-^ — 7T, a distinct from u and v. (10) 

Here the indices p and v refer to two distinct principal unit vectors and ej*, corresponding to the 
propagation and polarization direction, respectively. Clearly, these speeds are strictly within the range of 
microstability and causality, obeying 

< v^^, < 1. (11) 

Moreover, we prove in the appendix that the speeds of these modes are less than the speed of light also in 
generic propagation directions. 

Although the SUREOS is neither of the quasi-Hookean type discussed in paper I, nor of the closely 
related type introduced by Carter and Quintana|2|, we can still define an unsheared energy density p and 
shear modulus p, that are functions of the particle density n = nin2n^ only. This should in fact be possible 
for any physically reasonable equation of state of the type The procedure is simply to consider the 
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behaviour of the equation of state for smaU shear, i.e. when all linear particle densities are close to the 
unsheared value n^^^. To this end it turns out convenient to make the Misner type reparametrisation 



ni 



n2 = n^l'^e~^^+^^^ 



«3 



^1/3^-2:2-^3 a:i 



Consider first the quasi-Hookean equation of state of paper I for which the energy density is given by 

Pq.H. = P + 

where p and p. are functions of n whereas is given by 

1 



,s2 = 



12 

1 r 
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\n2 




V"3 


n2 ) 



"3 



smk? {\pi x\ — 3x2) + sinh^(2\/3 a^i) + sinh^(\/3a;i - 

Inserting the reparametrisation (|12|l into the equation of state Q instead gives 

p = An4/3[e-4-2 + g2(x2-y3.i) ^ g2(.2+y3.i)] + 5, 



3x2) 



(12) 

(13) 

(14) 
(15) 

(16) 



At fixed particle density n this expression for p can be identified with a cyclic three-particle Toda lattice 
potential reduced to the centre of mass system (c/. PH])- It directly follows that the zero shear state 
= 2:2 = is the only state that extremizes p (and hence also the energy per particle e = p/n) at fixed 
n. Moreover this extremum is a minimum if A is positive, which shall thus be assumed henceforth. Taylor 
expanding the two energy densities (|13|l and H16|) in x\ and X2 gives 



Pq.H. = P + 'Opix^ + Xi) + 0{\xt) 

p ^ 'iAn^l^^ + B + 12An^/^{x^ + x^) + 0{\x\^). 



(17) 
(18) 



Hence we can identify the SUREOS as a quasi-Hookean equation of state up to second order in x\ and X2, 
with unsheared energy density and shear modulus given by 



p = -iAn^l"^ + B 
p = 2An^l\ 



(19) 
(20) 



Clearly the shear modulus is very large, of the same order as the unsheared energy density unless B ^ iAn^l^ . 
As a consistency check, let us now use eqs. (|19|l and (|20|l to calculate the longitudinal and transversal wave 
speeds U|| and uj_ with respect to an isotropic background. To do so we first calculate the unsheared pressure 
and the bulk modulus: 



V 



dp 

dn 
dp 

dn 



P - 



An^/^ - B = \ip - AB) 



lAn^/\ 



Using the formulae of Carter [TT] (rederived in paper I), we find that the isotropic wave speeds are 



+ 

P + P 
P + P 



dp 
dp 



3p + p 



(21) 
(22) 

(23) 
(24) 



This is indeed consistent with our previous results, namely that the longitudinal wave speed is always equal 
to the speed of light while the principal transversal wave speeds are given by eq. H10|l . 
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3 Stationary solutions with rigid motion 



We here restrict our attention to stationary spacetimes, i.e. spacetimes admitting a timelike Killing vector 
We also assume that the elastic matter source of Einstein's equations undergoes rigid motion, by which 
we mean that the four- velocity m° is aligned with 1°" , i.e. 

= (25) 

for some scalar /. Since for elastic matter it always holds that 

ha''u''Tbc^ ha''u''{pUbUc+Pbc) ^0, hab ^ UaUb + gab, (26) 

it follows from Einstein's equations [T^ that the twist of <° is locally the gradient of some scalar oj, 

tOa = -tabcdt^^'t'^ = VaO^. (27) 

The four-dimensional physical spacetime (M, gab) is in this case completely determined by the three-dimen- 
sional structure (S, hab, f, i-^), where S is the three-manifold of Killing orbits and where hab, f and lu are fixed 
fields on E. By "fixed" we mean that these fields are all anihilated by £t, i-S- constant along the Killing 
orbits. Letting Da denote the connection associated with hab, the Einstein equations with stress-energy 
tensor Tab = pUaUb + Pab can be written in three-dimensional form as[12| 

np = _^ jr^D'^Lu DaLo (28) 

KPab = '^^^Gab + r'[iD'DJ)hab - DaDbf] + ^/^^ {D^UjD ,Uj)hab - 2DaUjDbUj]. (29) 

Since the energy density p and the three principal pressures (the eigenvalues of p'^b) all are constant along 
the Killing flowlines there will exist at least one functional relation, 

F(p,pi,p2,P3) =0 (30) 

for any stationary spacetime of this type. It is not difhcult to show that if the relation H30|l can be written 
such that it is invariant under permutations of the p^'s, then the spacetime represented by hab, f and to on 
S can, under quite general circumstances, be interpreted as a rigid motion solution to the Einstein equation 
for at least some elastic equation of state p — /9(ni, 712,713) and an associated material space metric kAB- 
Explicitly, if the relation = dp/ dn^ — p is inserted into eq. H3()|l . the latter becomes a partial differential 
equation for the dependent variable p(ni, 712,^3). If the mentioned permutation symmetry holds there will 
generically exist"*^, at least locally, a nonempty family of solutions /9(ni, 7i2, 713) that are invariant under 
permutations of the n^/s and hence can be interpreted as elastic equations of state of the type we focus 
on. Picking one member of the family, each linear particle density becomes an implicit function of three 
independent combinations of p and the p^'s, for instance by inverting the three relations = p^{ni, 712, 773). 
This in turn determines the tensor kab — X^j'i^i "^^^tiaRiib which by this construction trivially satisfies the 
necessary and sufficient conditions u'^kab — and Cukab = for it to be the puUback of a material space 
three-metric kAB from a three-manifold that can be identified with S. Since kAB should be positive definite, 
one must finally require that 77^ > for all p. It should be stressed that this method is in general of little 
practical use when it comes to finding exact solutions, since the equation of state and the material space 
metric will only be implicitly defined and cannot, except in special cases, be expressed in closed analytic form. 
However, one such special case is provided by the SUREOS, i.e. by eq. 0, for which eq. (|5n|) corresponds 
to the linear relation 0. That linear relation is identically satisfied for the SUREOS but in fact also for a 
wider class of equations of state that we do not discuss here. Taking thus the SUREOS to be the equation 
of state we can invert the relation to find that the linear particle densities are given by 



4 ^ {P ~ ~ 25) (p - - 25) 

2A{p-p^,-2B) ' ^ ' 



^We make no attempt to sort out under what precise conditions this will be true, since it is not our goal of this discussion. 
Such an analysis would probably be nontrivial since eq. I3UI can in principle be a complicated nonlinear function of its arguments. 
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where we have used eq. ® ■ Since p — — 2B are the eigenvalues of O^b = [p— 2B)h°-b — p"^b, the inverted 
formulae (|31|l implies that kab can be calculated from the stress-energy tensor via Q°'b according to 



[2^det(e^/)]-i/2 {e-^e,^ _ e-^e^, + ^[(e^,)^ - e^dO^dhab} , 



(32) 



where the last frame-invariant form, should it be cumbersome to diagonalize Pab, is more useful than the 
relation H31|) . 

The linear relation (PJ and Einstein's equations directly imply that the spacetime Ricci scalar R should 
take the constant value AkB. Since we must have ^ > for the equation of state to be physical, three 
inequalities are also imposed from 



p - Pf, - 2B = 2A{nf^+in 



fj,+2) 



> for all p., 



(33) 



which is the same as saying that Qab should be positive definite. All can subsequently be chosen to be 
positive, which is the last condition to be satisfied. Hence we have a very straightforward recipe for obtaining 
stationary rigid motion exact solutions to Einstein's equation with the SUREOS. We simply have to look 
for stationary spacetimes described by hab, f and lo as above, for which the four-dimensional Ricci scalar 



R = _ 21-^0" Da.f + ir^D^'ivDau; 



(34) 



is constant. This constant should then be identified with AkB and should hence be nonnegative if we want 
the equation of state parameter B to be. If, in addition, Qab is positive definite, we have a solution with 
respect to the material space metric (|32f) . The condition R = AnB may clearly be written as 



D^Daf 



2KB]f 



{r'^D'^uDaLO. 



(35) 



It would be very interesting to try to find twisting solutions, especially axisymmetric solutions that could 
describe rotating stellar models. However, we shall here focus on the simpler static case, Dai^ = 0, for which 
eq. (|35|l is linear and homogeneous in /. In this case eq. ()35(l is obviously particularly simple when the 
spatial Ricci scalar ^^^R is constant, which according to eq. H28|) corresponds to setting the energy density 
p constant. Morever, since ^^■'i? and p have the same sign we should have ^^^i? > in order to obtain a 
physically reasonable solution. The simplest case to consider is then to take hab to be locally isometric to 
the metric of a three-sphere of some radius rg, giving ^^^i? = 2k/3 — Q/r^ . Choosing coordinates such that 
hab takes the form 

dP = [dx^ + sin^x [dO^ + sin^^ dcj)"^)] (36) 

and expanding / according to 



smx 



1=0 m=-l 



we obtain the radial equation 



dx^ 



■ 2 

sm x 



Fi=0, 



(37) 



(38) 



where 

fc2 = 2(1 - Kr^B) (39) 

Whether the expansion (|37|l is appropriate to make depends on whether or not the resulting solution can 
be matched to an exterior asymptotically flat vacuum solution, or perhaps to some other matter solution 
of interest. This will be a nontrivial task to analyse, except for the spherically symmetric case; Qm — for 
I > 0, which we focus on in section 0] 
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3.1 Conformally rescaled formulation 



When dealing with a four-dimensional stationary metric it is often useful to work with the variables ^' = 
(minus the norm of the Killing vector t"), uj and the conformally rescaled spatial metric hab = ^hab, rather 
than /, UJ and hab- Making this reformulation here, we find that eq. (|35(l transforms into the equation 

which is completely written in terms of the metric hab] Da is its connection, R its Ricci scalar and indices are 
raised with its inverse, i.e. with h~^'^^ — ^"^/i*^^. Note that 5" and uj correspond to the real and imaginary 
part, respectively, of the Ernst potential £ = + iuj which in vacuum satisfies the equation 

We expect eq. to be a useful complement to eq. (^5]) . especially when looking for solutions that can be 
matched to exterior vacuum solutions, the latter being completely determined by giving the Ernst potential 
£ in the important case of axisymmetry. Again, a more detailed investigation of these matters is postponed 
for future work. 



4 Exact SSS solutions with constant energy density 

In the spherically symmetric case, the ansatz that the spatial metric hab be a three-sphere metric implies 
that the spacetime metric can be written as 

F \^ 2-F 



ds^ = dt^ + \dx^ + sin^x (dO^ + siii^e dcj}'^)] , = , (42) 

\sma; / ■' 2kB 

where F should satisfy eq. 1)38(1 with ^ = 0; 

F" = k^F. (43) 
We shall only discuss the solution that corresponds to a regular centre model, which reads 

= asinhfca;, (44) 

where a is a constant to be determined by the matching at the surface of the star. The energy density and 
principal (radial and tangential) pressures are given by 



P = p - Po (45) 

— (1 + cos^a;) sinh fcx -|- fc sin 2a; cosh fcx 

Pr = z—^ — — Po (46) 

6 sm X smh kx 

(cos^x + k^ sin^a;) sinh kx — ^k sin 2a; cosh kx 

Pt = ^ , . 2 . PO, (47) 

6 sm X smh kx 

whereas the dimensionless quotient m/r between the standard mass function m and the Schwarzschild radius 
r = ro sin x takes the simple form 

— i sin^x. (48) 
r ^ 

Since we are dealing with a regular centre solution, the radial and tangential pressure coincide at a; = 0. 
The value of this isotropic central pressure is given by 

Pe = i(2fc^-l)po = ^^-^S. (49) 
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No loss of generality is implied by assuming that k is non-negative (if real) so for positive B the physical 
range of k is 

-^<k<V2, (50) 
v2 

where the lower and upper limits correspond to zero and infinite central pressure, respectively. The k = V2 
solution is of course singular at the centre if S ^ 0, but regularity can be restored if B is set to zero, in 
which case the constant vq becomes a scaling parameter which is trivial in the sense that it has no relation 
to any scale in the equation of state. 

It is worth noting that it is possible to let k exceed \/2 if allowing B to take negative values, but this 
possibility will not be pursued here (but see fig. 

According to eq. the radial and tangential linear particle densities are given by 

iP-pt-2Br 

^ 2A{p-pr-2B) 



^ _ p-pr -2B 



(52) 



2A 

in terms of which the pull-back of the material space metric is 

kab = n^ran + n^tab, (53) 

where 

Ta = ro {dx)a, tab = siu^x {dO'^ + d4>'^)ab- (54) 
Obviously we could write down the metric H53|l explicitly in the coordinates (a;, 9, 0), but refrain from doing 
so since the result is not very illuminating. Note, however, that kab always has a regular centre whenever 
the spacetime metric gab does. 

From eqs. (|10|) . (|51|l and H52|) we can also calculate the three generally distinct transversal elastic wave 
speeds as 

2 



^ (55) 



t ^ "r 
2 TT-t 1 



(57) 



The surface of the model will be at a radius for which the radial pressure vanishes, allowing for a matching 
to the vacuum Schwarzschild solution to be made. For the allowed values of k such a surface always exists at 
a finite distance from the centre, although we have no closed analytic form for the solution to the equation 
Pr = 0. However, the constant a can be determined by the matching and is given by 

a = ^ sin 2a;s sinh fcxs, (58) 

where the subscript s refers to evaluation on the surface. 

The further properties of this constant energy density family of solutions are presented in figures ITlSl and 
their captions. 

4.1 Stability against radial perturbations 

According to the results of paper II, the radial perturbations for the constant energy density solutions are 
governed by the first order system 

= ro cos x\YC + P-^Tf\ (59) 

dx 

^ = -rQ cos x[{Q + uj''W )C + Yrf\, (60) 
dx 
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Figure 1: Mass-radius curve for the SUREOS constant energy density solution with 
B at a fixed positive value. The curve is parametrized by the dimensionless pa- 
rameter k in the range 1/V2 < k < \/2 or equivalently by the central pressure 
Pc = |i?(2fc^ — l)/(2 — A:^) taking all positive values. The direction of increasing k 
and Pc is the counterclockwise one. The maximum of the curve is at A; « 1.133 
and separates this one-parameter family of solutions into a stable branch with 
l/-\/2 < k < 1.133 and a branch with 1.133 < A: < having one unstable ra- 
dial mode. A peculiar feature of this mass-radius curve is that it not only starts out 
from the origin M = i? = at zero central pressure but also returns there in the 
limit of infinite central pressure. 



where the indendent variables ( and rj are closely related to the lagrangian perturbations of the Schwarzschild 
radius and the radial pressure, respectively. Moreover vq = — k'^)/{2KB) and 



2a^ (— cos 2x sinh kx + ^k sin 2x cosh kx) sinh^ A;x 



sm X cos X 



Y = - 



1 [(2 -I- fc^) sin^a; -|- 1] sinh kx — ^k sin 2x cosh kx 
ro (cos 2x sinh kx — i A; sin 2a; cosh A;a;) sin a; 



W = Pa- 



, tan^a; 
sinh^A;a; 



Q = Qi + Q2 

4 k'^ sin^x cosh^Ara; — cos 2a; sinh^A;a; 



Qi=P- 
Q2 = P{Z + Y^) 



sin^2a; sinh^A;a; 



(61) 

(62) 

(63) 
(64) 

(65) 
(66) 



Z = 



1 3 cos X [(2 cos^a; -|- 5) sinh^A;.T + A^ sin^x (5 cosh^ Ax — 3)] — A sin a: (1 -I- ^ cos^a; -|- A;^ sin^a;) sinh 2A;a; 



(cos 2a; sinh kx — hk sin 2x cosh A;a;) sin^a; cos x sinh Aa; 



(67) 
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A regular solution to the system (|59|l and (|60|) should leave the centre according to 



ri^CPo{3~Yo) + 0{x^), 



(68) 
(69) 



where C is an arbitrary constant and where Pq and Yq, as defined in paper II, are given by 



^0 



Yo 



3 kvq 
2. 



(70) 



(71) 



At the surface of the star the lagrangian perturbation of the radial pressure should vanish, implying 



We have numerically determined the first two frequencies ujq and uji of radial oscillations for B > and 
found that ujq changes sign from positive to negative at fc « 1.133, which to numerical precision^ coincides 
with the maximum mass model. As expected, is positive for all allowed values of k. The results are 
presented in figure |31 



^We use standard methods as described in 1131 . More precisely we find the frequencies using the shooting method employing 
a fourth-fifth order adaptive step size Runge-Kutta ODE-solver. The change of sign of the eigenvalue uJq coinsidos with the 
maximum of the mass to within more than six significant digits in the mass. 



= 0. 



(72) 
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0.5 




Figure 2: The compactness M/ R for the SUREOS constant energy solutfon, pfotted 
as a function of the dimensionless parameter k. Although the discussion in this paper 
is restricted to non- negative values of B, which requires that k < v^, the plot has 
been continued into the negative B sector k > \pl to display that the compactness 
then exceeds the perfect fluid Buchdahl limit value MjR = 4/9 for sufficiently large 
k. In fact, as k is taken to infinity the curve asymptotically approaches the horizon 
value MjR — 1/2. Such extreme models having a surface arbitrarily close to the 
horizon have previously been investigated for other matter types, for charged perfect 
fluids by de Fefice et al.^j EI- In the positive B sector l/-\/2 < k < y/2 the 
solutions are stable up to fc « 1.133, implying that the limiting compactness for the 
stable branch is M/R « 0.3408. Notably this value is slightly larger than 1/3 and 
hence there are stable solutions that are ultracompact in the sense that their surfaces 
are inside R = 3M at which the exterior Schwarzschild solution has closed lightlike 
geodesies. At fc = \/2 (infinite central pressure when B ^ 0) the compactness is 
M/R « 0.4096. 
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Figure 3: Energy density and pressure profiles for the marginally stable model with 
k w 1.133. The radius of this model is i? « 0.09857 x B'^/"^. The radial (tangential) 
pressure is monotonically decreasing (increasing) for other values of k as well. 
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Figure 4: Profiles of the speeds of clastic wave propagation for the marginally stable 
model with k « 1.133. While the single longitudinal wave speed v\\ is always unity 
(speed of light) the principal transversal wave speeds Vr±t, Vt±r and Vt±t start out from 
the unsheared value at the centre and diverges outwards as the pressure anisotropy 
increases. For smaller (larger) values of k, and thus of the central pressure, the profiles 
are qualitatively similar but the divergence is weaker (stronger). The transversal speed 
Vt±r is however always exactly equal to \/2. 
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The first two frequencies of radial oscillations 




1 .2 1 .3 1 .4 



Figure 5: The first two frequencies of radial oscillations for the family of constant 
energy density models with B > 0. The squared frequencies as well as the stellar mass 
M are plotted as functions of the dimensionless parameter k, displaying clearly that 
ujq tTirns negative precisely at the maximum mass model having k w 1.133. Although 
partially cut out from the plot, the square of the first overtone frequency u>i remains 
positive up to the infinite central pressure limit k = \pi. 
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5 A solution with nonconstant energy density 



The constant energy density of the family of solutions presented in the previous section conies about because 
the unsheared, i.e. compressional, energy density p decreases outwards at exactly the same rate as the 
shearing energy density a = p — p inreases. Here we will briefly present a different SUREOS solution which 
is more realistic in the sense that the energy density as well as both principal pressures are monotonically 
outwards decreasing, just to illustrate that such solutions do exist. The solution was found for B = only 
and although we have made some attempts to generalise it to nonzero B we have not been able to do so. 
In Schwarzschild coordinates the spacetime metric for the solution has the simple analytic form 

ds^ = sinc~^(r/ro) {-a^dt^ + dr^) + r'^d^f, sincx = (73) 

X 

To match the solution to a Schwarzschild exterior without rescaling the time coordinate, the constant a 
should be chosen as 

a = sinc^(i?/ro), (74) 

where R is the radius of the star, i.e. r = i? is the surface of vanishing radial pressure, to be determined 
below. Regularity at the centre is clearly guaranteed by the fact that sine a; = 1 at a; = 0. The energy 
density and principal pressures are given by 

K.r^ p = 1 + sinc^(r/ro) — 2sinc(r/ro) cos(r/ro) (75) 

nr^Pr = — 1 + 3 sinc^(r/ro) — 2 sinc(r/ro) cos(r/ro) (76) 
2m 

nr'^pt^ = 1 - sinc^(r/ro) (77) 

r 

(78) 



Pc = 5Pc = 2 • C^^) 

The profiles of these quantities are plotted in figureEl The radial pressure drops to zero very close to r = 2ro, 
or more precisely 

i?wl.9979ro, (80) 

implying a « 0.2075 and R/M « 2.5237. 

A similar analysis to that of section ETl reveals that the fundamental radial mode is unstable with squared 
frequency ujq « —0.07017 Tq"^ which for a typical neutron star with a radius of the order of 10 km corresponds 
to an e-folding time of about 0.4 ms. Clearly then, this solution is not very interesting on its own right, but 
if it could be generalised to a family with B being a free parameter we expect there to be a stable branch of 
solutions. 



6 Discussion and outlook 

Exact solutions to the Einstein equations can very rarely be used to completely replace numerically obtained 
solutions, with a few important exceptions such as the Kerr-Newman family of black hole solutions. The 
reason is of course that it is almost always impossible to find an exact solution that obeys all restrictions that 
are imposed by requiring that the solution be astrophysically realistic. However, if an exact solution satisfies 
at least some basic criteria for physicality, that solution has the potential of being a useful tool for gaining 
qualitative insights. Furthermore it can be very useful as a test background model when doing numerical 
studies of various types of perturbations. The constant energy density family of exact SSS solutions with 
stiff ultrarigid equation of state does, in our opinion, satisfy the relevant list of basic criteria that renders 
it useful in these respects. In fact we have already used it in the present paper to test that our analysis 
and numerical codes for radial oscillations of spherically symmetric elastic matter models are behaving as 
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0.05 0.1 0.15 0.2 

Figure 6: Energy density and pressure profiles for the B = solution with outwards 
decreasing energy density and eigenpressures. 

expected; Clearly it cannot be a coincidence that we find the zeroth mode turning unstable for the maximum 
mass member of the family, to numerical precision. So far this had only been tested for the numerically 
integrated models presented in paper I, which all have moderate pressure anisotropics, but now we know 
that it also works as expected for models with more extreme anisotropics. 

Since the presented method of generating SUREOS solutions does not rely upon spherical symmetry 
but only on stationarity and rigid motion, it would be very nice if one could also apply it to finding a 
rigidly rotating solution family. Although the equation of state is rather extreme, it is nevertheless always 
microstable and causal and hence, in our view, more physical than the equation of state p + 3p = po for 
Wahlquist's rotating perfect fluid solution^^l- The latter equation of state fails to be microstable since the 
squared speed of sound dp/ dp is obviously negative. 

Appendix 

Since the principal longitudinal speed is equal to the speed of light, it is of interest to find out whether 
this speed can become superluminal when the propagation vector z/° is perturbed around the principal unit 
vector e^. For the SUREOS the Fresnel tensor can be written in the simple form 

Qab ^ pj^ab ^ ^ab _ ^ab ^ ^g^) 

where L"'' is the tensor 

3 3 

(7=1 <y=l 
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Thus the characteristic equation, which given a unit propagation vector i'"' = '^a^a determines three 

propagation speeds v and polarization vectors can in this case be rewritten in the convenient form 

+ ^0. (83) 

Since aU three vectors ZJ^ are orthogonal to i^", so is the tensor L"'' which means that the characteristic 
equation always has the solution 

= 1 (84) 

l-a = l^a (85) 

Since this result holds regardless of the choice of propagation vector z/", we have proved the claim that 
the equation of state allows for a purely longitudinally polarized wave mode in all directions, i.e. not only 
principal directions which is the case generically. Moreover, the speed of these waves is direction independent 
and equal to the speed of light. Let us now consider L"'' as a tensor on the two-space orthogonal to (and, 
of course, to u°'). As such, its determinant (or, rather, the determinant of the mixed tensor L°'b) can be 
calculated according to 

\L'^^ L'^'^ ^ac^bd, ^ab — ^abc^'^ ■ (86) 

The result is 

3 3 

iL-'L-^eacebd ^Y.^P~' P.+i - '2B){p - p„+2 - 2B)u;^ = AA^n^ ^ n>,2. (87) 

cr=l (7=1 

Clearly, this scalar can only be zero if at least one linear particle density vanishes which only happens if the 
material space mapping is degenerate which is disallowed. Now, with being orthogonal to v'^ and of rank 
two, it follows that the only solution to the characteristic equation l|83|l with = 1 is the longitudinally 
polarized solution. For, if ia was to have a nonvanishing projection orthogonal to i/", that projection would 
have to be annihilated by L"'' which is impossible by the nondegeneracy property of the latter. As a corollary 
we thus have proved that the remaining two polarization modes always have propagation speeds less than 
unity. This follows from the facts that they should vary continuously with the propagation direction and 
that they are strictly less than unity for propagation in principal directions, in which case the polarization 
is purely transversal. 
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